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Observables of Lattice Gauge Theory in Minkowski Space 
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U(l) gauge fields are decomposed into a monopole and photon part across the phase transition from the 
confinement to the Coulomb phase. We analyze the leading Lyapunov exponents of such gauge field configurations 
on the lattice which are initialized by quantum Monte Carlo simulations. We observe that the monopole field 
carries the same Lyapunov exponent as the original U(l) field. Evidence is found that monopole fields stay 
chaotic in the continuum whereas the photon fields are regular. First results are presented for the full spectrum 
of Lyapunov exponents of the U(l) gauge system. 



1. Monopole and photon part of U(l) the- 
ory 

We are interested in the relationship between 
monopoles and classical chaos across the phase 
transition at j3 c w 1.01. We begin with a Ad U(l) 
gauge theory described by the Euclidean action 
S{Ui} = E P (1 " cos0 p ), where U t = U Xjfl = 
exp(i6 X:fl ) and 9 P = 9 X ^ + 9 X+ ^ V — 9 x+ p^ — 9 X>V . 
Following Ref. [1], we have factorized our gauge 
configurations into monopole and photon fields. 
The U(l) plaquette angles 9 X41V are decomposed 
into the "physical" electromagnetic flux through 
the plaquette 9 XifJjL , and a number m Xjfil/ of Dirac 
strings through the plaquette 



ing zero values of these exponents. Finally neg- 
ative exponents belong to irrelevant directions 
in the phase space: perturbation components 
in these directions die out exponentially. Pure 
gauge fields on the lattice show a characteris- 
tic Lyapunov spectrum consisting of one third of 
each kind of exponents [2]. Assuming this gen- 
eral structure of the Lyapunov spectrum we first 
investigate its magnitude, namely the maximal 
value of the Lyapunov exponent, L max . 

The general definition of the Lyapunov expo- 
nent is based on a distance measure d(t) in phase 
space, 



2-tt m 1)P 
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L := lim lim —In -. \ . 

t^oo<j(o)^ot d(Q) 



(2) 



where 9 X ^ V G {—it, +7t] and m 1#l/ ^ is called a 
Dirac plaquette. 

2. Classical chaotic dynamics from quan- 
tum Monte Carlo initial states 

Chaotic dynamics in general is characterized 
by the spectrum of Lyapunov exponents. These 
exponents, if they are positive, reflect an expo- 
nential divergence of initially adjacent configu- 
rations. In case of symmetries inherent in the 
Hamiltonian of the system there are correspond- 



In case of conservative dynamics the sum of all 
Lyapunov exponents is zero according to Liou- 
ville's theorem, Li = 0. We utilize the gauge 
invariant distance measure consisting of the lo- 
cal differences of energy densities between two 3d 
field configurations on the lattice: 



d := — V \trU P - trU'p 
N P 1 F 



(3) 



Here the symbol J^p stands for the sum over all 
Np plaquettes, so this distance is bound in the 
interval (0, 2 AO for the group SU(N). U P and U' P 
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Figure 1. Exponentially diverging distance in real time of initially adjacent U(l) field configurations on 
a 12 3 lattice prepared at = 0.9 in the confinement (left) and at — 1.1 in the Coulomb phase (right). 



are the familiar plaquette variables, constructed 
from the basic link variables U x i, 



U x .i = exp (aA x ^T) 



(4) 



located on lattice links pointing from the posi- 
tion x = (xi,X2,x 3 ) to x + aei. The generator 
of the group (7(1) is T = —ig and A x ^ is the 
vector potential. The elementary plaquette vari- 
able is constructed for a plaquette with a cor- 
ner at x and lying in the ij-plane as U x ,ij = 
Ux.iUx+ijUl^^l j. It is related to the mag- 
netic field strength B x y. 



(5) 



U x>i j = exp ( , I Irj. l ) 

The electric field strength E x ^ is related to the 
canonically conjugate momentum P x .i — aU x ,i 
via 
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The Hamiltonian of the lattice gauge field sys- 
tem can be casted into the form 



H 



E 



(p,p) 
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(U,V) 



(7) 



Here the scalar product stands for (A, B) = 
Re(AE?t). The staple variable V is a sum of triple 
products of elementary link variables closing a 
plaquette with the chosen link U . This way the 
Hamiltonian is formally written as a sum over link 
contributions and V plays the role of the classical 
force acting on the link variable U. 



We prepare the initial field configurations from 
a standard four dimensional Euclidean Monte 
Carlo program on a 12 3 x 4 lattice varying the 
inverse gauge coupling [3] . We relate such four 
dimensional Euclidean lattice field configurations 
to Minkowskian momenta and fields for the three 
dimensional Hamiltonian simulation by selecting 
a fixed time slice of the four dimensional lattice. 

3. Chaos, confinement and continuum 
limit 

We start the presentation of our results with 
a characteristic example of the time evolution of 
the distance between initially adjacent configura- 
tions. An initial state prepared by a standard four 
dimensional Monte Carlo simulation is evolved 
according to the classical Hamiltonian dynamics 
in real time. Afterwards this initial state is ro- 
tated locally by group elements which are chosen 
randomly near to the unity. The time evolution of 
this slightly rotated configuration is then pursued 
and finally the distance between these two evolu- 
tions is calculated at the corresponding times. A 
typical exponential rise of this distance followed 
by a saturation can be inspected in Fig. 1 from 
an example of U(l) gauge theory in the confine- 
ment phase and in the Coulomb phase. While 
the saturation is an artifact of the compact dis- 
tance measure of the lattice, the exponential rise 
(the linear rise of the logarithm) can be used for 
the determination of the leading Lyapunov expo- 
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Figure 2. Lyapunov exponent of the decomposed 
U(l) fields as a function of coupling. 



nent. The left plot exhibits that in the confine- 
ment phase the original field and its monopole 
part have similar Lyapunov exponents whereas 
the photon part has a smaller L max . The right 
plot in the Coulomb phase suggests that all slopes 
and consequently the Lyapunov exponents of all 
fields decrease substantially. 

The main result of the present study is the de- 
pendence of the leading Lyapunov exponent L max 
on the inverse coupling strength (3, displayed in 
Fig. 2. As expected the strong coupling phase 
is more chaotic. The transition reflects the crit- 
ical coupling to the Coulomb phase. The plot 
shows that the monopole fields carry Lyapunov 
exponents of nearly the same size as the full U(l) 
fields. The photon fields yield a non-vanishing 
value in the confinement ascending toward (3 = 
for randomized fields which indicates that the de- 
composition works well only around the phase 
transition. 

An interesting result concerning the continuum 
limit can be viewed from Fig. 3 which shows the 
energy dependence of the Lyapunov exponents for 
the U(l) theory and its components. One ob- 
serves an approximately linear relation for the 
monopole part while a quadratic relation is sug- 
gested for the photon part in the weak coupling 
regime. From scaling arguments one expects a 
functional relationship between the Lyapunov ex- 
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Figure 3. Comparison of average maximal Lya- 
punov exponents as a function of the scaled av- 
erage energy per plaquette. The full U(l) theory 
shows an approximately quadratic behavior in the 
weak coupling regime which is more pronounced 
for the photon field. The monopole field indicates 
a linear relation. 



ponent and the energy [2,4] 
L(a) cx a k - 1 E k {a), 



(8) 



with the exponent k being crucial for the contin- 
uum limit of the classical field theory. A value of 
k < 1 leads to a divergent Lyapunov exponent, 
while k > 1 yields a vanishing L in the contin- 
uum. The case k = 1 is special allowing for a 
finite non-zero Lyapunov exponent. Our analy- 
sis of the scaling relation (8) gives evidence, that 
the classical compact U(l) lattice gauge theory 
and especially the photon field have k « 2 and 
with L(a) — ► a regular continuum theory. The 
monopole field signals k « 1 and stays chaotic 
approaching the continuum. 



4. Spectrum of the stability matrix 

Instead of the classical determination of the 
Lyapunov exponent, which has been calculated 
by the rescaling method in the preceding sections, 
we now use the monodromy matrix approach [5] . 
The Lyapunov spectrum Lj is expressed in terms 
of the eigenvalues Aj of the monodromy matrix 
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Figure 4. Full complex eigenvalue spectrum of the monodromy matrix for t — at (3 — 0.9 in the 
confinement (left) and at /3 = 1.1 in the Coulomb phase (right). 

M: 



1 /•' 

= lim - / A^t'Jdt 7 , *=1,...,/ , (9) 

where A, (i) are the solutions of the characteristic 
equation for / degrees of freedom 



det (Ai(t)l-M(t)) =0 



(10) 



at a given time t. Here M is the linear stability 
matrix, 



M = 



( dtr ou \ 

dU dP 



dP dP 
\ dU dP I 
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beta = 0.5 
beta = 0.9 



beta, 
beta. 



(11) 



with the three dimensional lattice action S3. 

Figure 4 displays the complex eigenvalues for 
selected configurations prepared by a quantum 
Monte Carlo heat-bath algorithm. In the con- 
finement phase ((3 = 0.9) the eigenvalues lie on 
either the real or on the imaginary axes. This is a 
nice illustration of a strange attractor of a chaotic 
system. Positive Lyapunov exponents eject the 
trajectories from oscillating orbits provided by 
the imaginary eigenvalues. Negative Lyapunov 
exponents attract the trajectories keeping them 
confined in the basin. In the Coulomb phase 
{fi = 1.1) the real Lyapunov exponents become 
rare to eventually vanish in the continuum limit 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
i/N LINK 



Figure 5. Real part of the squared spectrum of 
the stability matrix decreasing with increasing (3. 



(see Fig. 3). In all cases the spectrum is symmet- 
ric with respect to the real and imaginary axes: 
the former property is due to the fact that the 
equations of motion are real, the latter is due to 
the Hamiltonian being conservative (time inde- 
pendent) . Also a number of zero- frequency modes 
occur which is connected to symmetry trans- 
formations commuting (in the Poisson bracket 
sense) with the Hamiltonian, such as time inde- 
pendent gauge transformations. 
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Figure 6. Power spectrum at /3 = (left) and at /3 

Figure 5 shows the real part of the squared 
Lyapunov spectrum for several couplings. The 
overall pattern is similar to that obtained earlier 
for SU(2) studies with random configurations [5]. 
It exhibits an increasing number of zero modes 
with increasing (3. Figure 6 gives insight into the 
distribution of the squared eigenvalues which is 
symmetric for (3 = and is shifted to negative 
values towards the Coulomb phase. 

5. Summary 

This contribution was devoted to the relation 
between monopoles, Lyapunov exponents and 
confinement in a (classical) compact U(l) the- 
ory. A decomposition into a monopole and pho- 
ton part showed that the monopole field carries 
the same Lyapunov exponent as the original U(l) 
field. In the continuum limit the photon and U(l) 
fields behave regularly while the monopole part 
stays chaotic. This was substantiated by the spec- 
trum of the monodromy matrix exhibiting an in- 
terplay between positive, imaginary and negative 
Lyapunov exponents in the confinement phase 
changing to a pure imaginary spectrum deep in 
the Coulomb phase leading to the regularity of 
the Maxwell theory. It will be interesting to com- 
pute the spectrum of the stability matrix on the 
photon and monopole part of the U(l) field. An 



100 



80 




1.1 in the Coulomb phase (right). 

analysis of the behavior of the monopoles in the 
classical U(l) theory in Minkowski space is a chal- 
lenge by itself. 
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